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ABSTRACT EDITOR'S SUMMARY 

Larger brains tend to have more folded cortices, but what makes the cortex fold has remained unknown. We 
show that the degree of cortical folding scales uniformly across lissencephalic and gyrencephalic species, 
across individuals, and within individual cortices as a function of the product of cortical surface area and the 
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square root of cortical thickness. This relation is derived from the minimization of the effective free energy 
associated with cortical shape according to a simple physical model, based on known mechanisms of axonal 
elongation. This model also explains the scaling of the folding index of crumpled paper balls. We discuss the 
implications of this finding for the evolutionary and developmental origin of folding, including the newfound 
continuum between lissencephaly and gyrencephaly, and for pathologies such as human lissencephaly. 



The expansion of the cerebral cortex, the most obvious feature of 

Related Resources 

mammalian brain evolution, is generally accompanied by increasing 

degrees of folding of the cortical surface into sulci and gyri (1_). Cortical in Science Magazine 

folding has been considered a means of allowing numbers of neurons in PERSPECTIVE 
the cerebral cortex to expand beyond what would be possible in a brain evolution 

lissencephalic cortex, presumably as the cortical sheet expands laterally Knowing when to fold them 
with a constant number of neurons beneath the surface (2, 3). Although Georg F. Striedter, Shyam Srinivasan 
some models have shown cortical convolutions to form as a result of Science 3 July 2015: 31-32. 

cortical growth (4, 5), the mechanisms that drive gyrification remain to 
be determined, and the field still lacks a mechanistic and predictive, 
quantitative explanation for how the degree of cortical folding scales 

across species. Moreover, recent systematic analyses of cortical folding have made clear that gyrification actually 
scales differently across mammalian orders, across clades within an order, and across individuals as a function of 
increasing brain volume (6-9). These apparent discrepancies have led to the view that different mechanisms must 
regulate cortical folding at the evolutionary, species-specific, and ontogenetic levels (7). 

We undertook a systematic analysis of the variation in cortical folding across a large sample of mammalian 
species in search of a universal, unifying relationship between cortical folding and morphological properties of the 
cerebral cortex. We examined two data sets: our own, which includes numbers of cortical neurons and cortical 
surface areas (10-21). and another consisting of published data on cortical surface area, thickness, brain volume, 
and folding index, but not numbers of cortical neurons 22-24) (table S1). 

In the combined data set, there is a general correlation between total brain mass and the degree of cortical folding, 
and the two data sets overlap in their distribution ( Fig. 1A , compare black and colored data points). However, the 
power function that relates the folding index of gyrencephalic species to brain mass has a fairly low r 2 and a 95% 
confidence interval that excludes many species ( Fig. 1A ). Striking and well-known outliers in this relationship are 
cetaceans (as a whole) and the manatee, but the capybara, the greater kudu, and humans also lie outside of the 
confidence interval ( Fig. 1A ). This indicates that cortical folding is not a homogeneous function of brain mass. 
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Fig. 1. Scaling of cortical folding index and total cortical 
surface area. 

Data points in black are taken from the literature; points in colors 
are from our own data set, except for cetaceans. (A to E) Folding 
index scales across all gyrencephalic species in the combined 
data sets as power functions of (A) brain mass, with exponent 
0.221 ± 0.01 8 (r 2 = 0.751 , P < 0.0001 ); (B) number of cortical 
neurons, with exponent 0.1 68 ± 0.032 (r 2 = 0.573, P < 0.0001 ; 



Download PowerPoint Slide for Teaching j not plotted); (D) total cortical area, with exponent 0.257 ± 0.014 

1 (r 2 = 0.872, P< 0.0001); and (E) average cortical thickness, with 
a nonsignificant exponent (r 2 = 0.054, P = 0.1430; not plotted). 
(C) Total cortical surface area of the cerebral cortex scales across primate species with an exponent of 0.91 1 
± 0.083 (r 2 = 0.938, P < 0.0001 ) and across nonprimate species with an exponent of 1 .248 ± 0.037 (r 2 = 
0.989, P< 0.0001). (F) Total cortical surface area varies across lissencephalic species as a linear function of 
the exposed surface area, but as a power function with an exponent of 1 .242 ± 0.01 8 across noncetacean 
gyrencephalic species (r 2 = 0.992, P < 0.0001 ). Dashed lines are 95% confidence intervals for the fitted 
functions. Dashed lines are 95% confidence intervals for the fitted functions. 



Although all cortical hemispheres with fewer than 30 million neurons are lissencephalic in our data set, and the 
correlation between folding index and number of neurons is significant across gyrencephalic species (Spearman 
correlation, p = 0.7741 , P < 0.0001 ), the degree of gyrification is much larger in artiodactyls than in primates for 
similar numbers of cortical neurons ( Fig. 1B ). Additionally, the elephant cortex is about twice as folded as the 
human cortex, although the former has only about one-third the number of neurons found in the latter ( Fig. 1B , "e" 
and "h"). The cortical surface area across species expands sublinearly with the number of cortical neurons in 
primates and supralinearly in other species ( Fig. 1C ). As a consequence, the average number of neurons per mm 2 
of cortical surface is highly variable across species, ranging in our data set from 1 0,752 in the African elephant 
( 15 ) to 138,606 in the squirrel monkey (20). Cortical expansion and folding are therefore neither a direct 
consequence of increasing numbers of neurons nor a requirement for increasing numbers of neurons in the cortex. 

In comparison to the poor fit between folding index and total brain mass ( Fig. 1A ), a better fit is found for total 
surface area of the cerebral cortex in the two data sets ( Fig. 1D ). In this case, there is better overlap across 
afrotherians, glires, primates, and artiodactyls, although cetaceans, the manatee, and humans are still major 
outliers. Interestingly, all species with a cortical surface area below 400 mm 2 are lissencephalic in the two data 
sets. Similarly, all species with average cortical thickness below 1 .2 mm are lissencephalic, but the folding index 
does not vary as a significant power function of cortical thickness across gyrencephalic species ( Fig. 1E ). 



The folding index shows a sharp inflection between smooth and gyrated cortices, so it is unlikely that a universal 
model in terms of this variable alone could be derived. Because the folding index is the ratio of total surface area 
Aq to exposed surface area A^, we next examined directly how A^ scales with Aq ( Fig. 1F ). In the combined data 
sets, for the species with small Aq (<400 mm 2 ) there is no folding, such that A^_ equals Aq ( Fig. 1F , green line). 
This linear relationship extends to the manatee cerebral cortex, even though its Ag is much larger than 1000 mm 2 . 
In contrast, for all noncetacean gyrencephalic species, Aq increases with /\ E 1 -242±0.018 ^.2 _ 0.992, p < 0.0001), 
significantly above linearity ( Fig. 1F , red line), meaning that as total surface area increases, it becomes 
increasingly folded. Cetaceans fall above the 95% confidence interval of the function, which indicates that these 
cortices are more folded than similarly sized cortices in noncetaceans. 

The finding that Aq scales as a power law of Ag means that gyrification is a property of a cortical surface that is 
self-similar down to a fundamental scale (the limit area between lissencephaly and gyrencephaly). This strongly 
suggests the existence of a single universal mechanism responsible for cortical folding (the alternative being some 
improbable multiscale fine-tuning) that over a range of scales generates self-similar, or fractal, surfaces. 

Fractals can be characterized by the power-law scaling between intrinsic and extrinsic measures of an object's 
size, such as Aq and A^. In this case, the fractal dimension d of the cortical surface is twice the value of the 
exponent relating Aq to Ag (given that Ag in turn scales with the square of the linear dimension of the cortex). 
Given that Aq scales with / 4 E 1 -242±0.018 across noncetacean gyrencephalic brains, then d = 2.484 ± 0.036. This 
value is remarkably close to the fractal dimension 2.5 of crumpled sheets of paper (25). which are fractal-like self- 
avoiding surfaces thin enough to fold under external compression while maintaining structural integrity. 

Empirically, we conceive the fractal folding (or lack thereof) of the cortical surface as a consequence of the 
minimization of the effective free energy of a self-avoiding surface of average thickness T that bounds a volume 
composed of fibers connecting distal regions of said surface. Our model incorporates the known mechanics and 
organization of elongating axonal fibers (26, 27), as described in the supplementary materials. It predicts that from 
a purely physical perspective, Aq, Ag, and Tare related by the power law T^ I2 Aq = kA^ IA . (The exponent 5/4 is 
the only value for which the constant k is adimensional.) 

We first tested whether our model based on the minimization of the effective free energy of a self-avoiding surface 
could explain the well-known fractal folding of a self-avoiding surface: paper. We examined how the exposed 
surface area of crumpled paper balls, Ag, scales with increasing total surface area, Aj, and thickness, T, of office 
paper (in this case, under forces applied externally by the experimenter's hands). As shown in Fig. 2A , Aj ~ 
Ag 1 ■ 234±0 033 for crumpled single sheets, a value similar to that for gyrencephalic cortices. Increasing T (by 
stacking sheets before crumpling) displaces the curves to the right ( Fig. 2A ) but leaves their slope largely 
unaltered, resulting in similar-looking but less folded paper balls ( Fig. 2B ). However, the product T 1/2 Aj varies 
proportionately to /\ E 1 - 105±0 °22 as a single, universal power function across all paper balls of different surface 



areas and thicknesses ( Fig. 2C ), as predicted by our model. This conformity indicates that the coarse-grained 
folding of a sheet of paper subjected to external compression depends simply on a combination of its surface area 
and thickness. 



Fig. 2. The degree of folding of crumpled paper balls is a 
function of surface area and thickness as predicted by our 
model. 

(A) Relationship between total surface area of A4 to A1 1 sheets 
of office paper and the exposed surface area of the crumpled 
sheet of paper, with a power function of exponent 1 .234 + 0.033 
for a single sheet of thickness 0.1 mm. (B) Increasing the 
thickness of the paper to be crumpled by stacking two to eight 
sheets displaces the curves to the right, that is, decreases the 
folding index of the resulting paper balls. (C) However, all 
crumpled paper balls of varying total surface area and thickness 
exhibit the same relationship, with the product T^ l2 Aj varying 
proportionately to -10510.022 (r 2 = 0 g83 p < 0 0001 y Co]or 

gradations correspond to thickness in millimeters, as shown in 
each panel. 



We next examined whether our model predicts the folding of the mammalian cerebral cortex by plotting the product 
T 1/2 /Ag as a function of for the combined data sets. This yielded a power function with an exponent of 1 .329 ± 
0.014, with a very high r 2 of 0.996 for the noncetacean gyrencephalic species in the combined data sets ( Fig. 3A , 
red line). Note that this function, although calculated for gyrencephalic species, overlaps with lissencephalic 
species. Including lissencephalic species (but still excluding cetaceans) actually improved the fit, with r 2 = 0.998, 
and yielded an exponent of 1 .305 ± 0.01 0, which is close to the expected value of 1 .25. Adding cetaceans to the 
analysis resulted in a small change of the fit ( Fig. 3A , black line). Remarkably, the function fitted exclusively for 
lissencephalic species also predicted the relationship between T^ I2 Aq and A^ in gyrencephalic species ( Fig. 3A , 
green line) — and species such as the manatee and other afrotherians are no longer outliers. 

Fig. 3. The degree of folding of the mammalian cerebral cortex is a single function of surface area and 
thickness across lissencephalic and gyrencephalic species alike, although thickness scales as order- 
specific functions of cortical surface area. 

(A) The product T 1/2 /A G varies with /A E 1 ■ 32Q±00U (i 2 = 0.996, P < 0.0001 ) across noncetacean gyrencephalic 
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species in the combined data set (red line), with /\ E 1 - 32 5±0.009 ^2 
= 0.997, P< 0.0001, k = 0.157 ± 0.012) across all species 
(including cetaceans; black line), and with /\ E 1 -292±0.027 ^2 _ 
0.994, P < 0.0001) across lissencephalic species alone (green 
line). Note that the function plotted for lissencephalic species 
predicts the product T^ I2 Aq for gyrencephalic species equally 
well as the functions plotted for gyrencephalic species 
themselves. (B) Cortical thickness varies with cortical surface 
area ^0.555+0.053 ^ _ 0 .887, P < 0.0001 ) across lissencephalic 
species in the combined data set (green line), but with 
^ G 0.160±0.025 (f 2 = 0.703, P< 0.0001) across primates (red line), 
and with a q 0.334±0.072 (r 2 = 0 879 p = 0 .0185) across artiodactyl 
species (pink line). All fits exclude cetaceans. Dashed lines 
I indicate the 95% confidence intervals for the fitted functions. 



Given the theoretical relation T^ I2 Aq = kA^ 25 , it follows that lissencephalic species (for which Aq = A^) are those 
that meet the condition T= k 2 A(^ t2 . In contrast, all species for which T< k 2 A(± 12 are predicted to be 
gyrencephalic, with Aq > A^ (the alternative where T> k 2 A(^ 12 would result in Aq < A^, which is geometrically 
impossible). Indeed, in the combined data set, we find that ja Aq°- 555±0 053 (p< 0.0001) across lissencephalic 
species ( Fig. 3B , green line). All gyrencephalic species data points fall to the right of the lissencephalic 
distribution; that is, their Aq values are larger than predicted for a cortical thickness that would allow lissencephaly. 
The precise relationship between T and Aq across gyrencephalic species differs across orders, with a much 
smaller exponent for primates than for artiodactyls ( Fig. 3B , red and pink lines). Thus, within the single universal 
relationship that describes cortical expansion, there is a transition point between smooth and folded cortices: 
Gyrencephaly ensues when Aq expands in area faster than T 2 . For gyrencephalic species, the rate of expansion 
of cortical thickness relative to expansion of the cortical surface varies across orders, but the product P I2 Aq still 
varies as a universal function of/A^ 1,26 ioA^ 33 . 

We also found the same universality between the product T^ I2 Aq and A^ across coronal sections along the 
anteroposterior axis of the cortical hemisphere of a single individual, of different individuals, and even different 
species ranging from small rodents to human and elephant (fig. S1). 



The finding that A^ scales across all lissencephalic and gyrencephalic mammals (and even across species usually 
regarded as outliers such as the manatee and cetaceans) as a single power law of 1 V2 Aq indicates that 
gyrification is an intrinsic property of any mammalian cortex. Further, because the degree of folding can be 



described by the simple equation generated by our model (which also applies to crumpled sheets of paper), folding 
must occur as it minimizes the effective free energy of the cortical surface. Folding is therefore an intrinsic, fractal 
property of a self-avoiding surface, whether biological or not, subjected to crumpling forces. As such, this scaling of 
cortical folding does not depend on numbers of neurons or how they are distributed in the cortical sheet, but simply 
on the relative lateral expansion of this sheet relative to its thickness, regardless of how densely neurons are 
distributed within it. 

The finding that cortical folding scales universally across clades, species, individuals, and parts of the same cortex 
implies that the single mechanism based on the physics of minimization of effective free energy of a growing 
surface subject to inhomogeneous bulk stresses applies across cortical development and evolution. This is in stark 
contrast to previous conclusions that different mechanisms regulated folding at different levels (7); such 
conclusions may reflect the traditional emphasis on the relationship between folding degree and brain volume 
8), which is indeed diverse across orders, across species, and across individuals of a same species (6, 8). Also, 
the dependence of cortical folding on a simple combination of Aq and T implies that any alterations, such as 
defects in cell migration, that lead to increased Tor decreased Aq (or both) are expected to decrease cortical 
folding, exactly as found in human pathological lissencephaly (28). This might also be the case for the 
lissencephalic brain of birds, where a very thick telencephalon of small surface area surrounds the subpallial 
structures. 

Finally, our findings indicate that cortical folding did not evolve, in the sense of a new property specific to some 
clades but not others. Similarly, there is no such thing as "secondary lissencephaly" (29), nor are there two clusters 
of gyrencephaly (9). Rather, what has evolved, we propose, is a faster increase in Aq relative to T 2 in development 
— and at different rates in different mammalian clades, which thus become gyrencephalic at different functions of 
Aq or different numbers of neurons. 

Remarkably, there is no a priori reason for lissencephaly, considering that/Ag and T ultimately result from different 
biological processes: lateral expansion of the progenitor cell population for Aq, radial neurogenesis and cell growth 
for TOO). Similarly, there is no a priori reason for the cortex to become gyrencephalic once past a certain surface 
area — unless the rate of (lateral) progenitor cell expansion inevitably outpaces the rate of (radial) neurogenesis at 
this point, which apparently occurs typically when Aq reaches 400 mm 2 . We propose that, starting from the earliest 
and smallest (and smooth) mammalian brains (31), the cortical surface initially scaled isometrically, with Aq ~ J 2 . 
Gyrencephaly ensued in each clade as soon as this lockstep growth changed, with Aq now increasing faster than 
T 2 . Probable mechanisms involved are those that control the rate of neurogenesis and increases in cell size 
relative to the rate of progenitor and intermediate progenitor cell proliferation in early cortical development. Rapid 
increases in numbers of intermediate progenitor cells would lead to gyrencephaly, although not through the 
generation of larger numbers of neurons, as previously thought (7, 30, 32), but rather through the simple lateral 
expansion of the resulting cortical surface area at a rate faster than the cortical thickness squared. 
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